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Abstract: Phylogenetic relationships of Arundinaria and related genera (Pleioblastus, Pseudosasa, Oligostachyum, Bashania, Clavinodum, 
etc.) were assessed by analyzing the sequences of the nrDNA internal transcribed spacer (ITS) and the cpDNA trnL-F intergenic spacer 
(IGS). Comparison with trnL-F IGS sequence, the ITS region provided the higher number of parsimony informative characters, and the 
interspecific variation of the ITS sequence was higher than that of the trnL-F IGS sequence.The tree obtained by combining both sets of data 
showed that the species sampled in Arundinaria and the related genera were monophyletic and divided into two clades. The relationships 
and positioning of all the taxa surveryed (including A. oleosa, A. hsienchuensis, A. chino, A. amara, A. yixingensis, A. amabilis, A. fortunei, 
A. pygmaea, A. gramineus, A. fargesii, A. faberi, A. hupehense, Pseudosasa japonica cv. Tsutsumiana, P. japonica and Brachystachyum 
densiflorum ) were also discussed. The results from the sequences were broadly consistent with morphological characters, appearing all these 
taxa sampled belong to the genus of Arundinaria. The topologies of the trees generated from individual data and the combined data were 
similar. 
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Introduction 

Arundinaria is the earliest genus to be established in the 
genera of monopodial bamboos with the most diverse bamboos 
species in the world. It is reported that there were about 400 
species under Arundinaria. However, up to now, there is no 
comprehensive system of classification which indicates species 
relationships within the bamboos in the Arundinaria and the 
related genera (including Pleioblastus, Pseudosasa, 
Oligostachyum, Bashania, Clavinodum, etc.) due to systematic 
difficulties of the availability of good taxonomic characters (Guo 
& Li 2002). The systematical classification of Arundinaria is still 
controversial, with taxa which have been alternatively merged 
into the genus of Arundinaria, or segregated into several related 
genera (e.g. Pleioblastus, Pseudosasa, Oligostachyum, Bashania, 
Clavinodum, etc.). Different views on the status of these genera 
have been made by Yang et al. (1993), Soderstrom (1987), Li et 
al. (2002), etc. 

Molecular data from DNA studies have been widely applied in 
genetic variation, systematical classification, and phylogenetic 
relationship of bamboos, which mainly based on randomly 
amplified polymorphic DNA (RAPD) (Ding 1998), simple 
sequence repeat (SSR) (Li et al. 2002), and restriction fragment 
length polymorphism (RFLP) (Friar & Kochert 1994), etc. The 
studies based on two different regions of the bamboo genome, 
ITS sequences from nrDNA and trnL-F sequences from cpDNA, 
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have been scarcely reported before. The objectives of the this 
study are: (1) to evalute the level at which molecular data from 
these two regions provide phylogenetic information, and (2) to 
test if the combination of the ITS region and the trnL-F 
intergenic spacer (IGS) sequences could provides valuable data 
for phylogenetic relationship within the species traditionally 
described to the Arundinaria and the related genera (including 
Pleioblastus, Pseudosasa, Oligostachyum, Bashania, 
Clavinodum, etc.). 

Materials and methods 

Seventeen species representing related genera of Arundinaria 
(Arundinaria gramineus is the model species of Pleioblastus, 
Arundinaria fargesii the model species of Bashania, Arundinaria 
sulcata the model species of Oligostachyum, Pseudosasa 
japonica the model species of Pseudosasa, Arundinaria 
oedogonata is the model species of Clavinodum ) and one 
outgroup ( Phyllostachys edulis) were sampled (Table 1). 

Total DNA was extracted from fresh leaves by the DNeasy 
Plant Mini Kit (QIAGEN). Nuclear ribosomal internal transcribed 
spacers (ITSL 5.8S, ITS2 region ) were amplified using the forward 
primer (ITS-1 at 18S: 

5'-AGAAGTCGTAACAAGGTTTCCGTAGG-3') and the 
reverse primer (ITS-4 at 26S: 

5 ’-TCCTCCGCTTATTG ATATGC-3'), the trnL-F intergenic 
spacer of the chloroplast genome was amplified using the primers 
e and f (Taberlet et al. 1991).The PCR products were purified by 
Wizard PCR Preps DNA Purification System (Promega), and 
then sequenced directly on AB1 310 DNA automated sequencer 
using ABI Prism Bigdye™ Terminator Cycle Sequencing Beady 
Reaction Kit. 

The entire ITS and trnL-F sequences of 18 taxa were aligned 
with the CLUSTAL W multiple sequence alignment program 
(Thomson et al. 1997). and were adjusted manually in the 
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sequences where necessary. Phylogenetic relationships were 
analysed by PAUP 4.0 (Swofford 1998) with all changes 
weighted equally, using heuristic search, random addition 
sequence, and tree bisection-reconnection swapping. Gaps were 


treated as missing data. Searches were conducted on the ITS data, 
on the trnL-F data and on the ITS and the trnL-F combined data. 
Bootstrap analyses of 1 000 replications were performed to show 
relative supports for individual clades. 


Table 1. Materials of 18 species of the bamboos sampled in this report 


Taxa 

Sources 

Collector 

Voucher® 

Arundinaria oleosa 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ2—103 

Arundinaria amara 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ15—03 

Arundinaria hsienchuensis 

Hangzhou Bamboo Garden 

ZHUGE Qiang 

HZS—12 

Arundinaria yixingensis 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ15—07 

Arundinaria gramineus 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ4—36 

Arundinaria chino 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ15—10 

Pseudosasa japonica 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ16—03 

Pseudosasa japonica cv. Tsutsumiana 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ9—15 

Arundinaria amabilis 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ6—15 

Arundinaria fortunei 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ5—05 

Arundinaria pygmaea 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ8—10 

Arundinaria sulcata 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ3—06 

Arundinaria hupehense 

Hangzhou Bamboo Garden 

ZHUGE Qiang 

HZ4a—01 

A rundinaria fargesii 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ3—20 

Arundinaria oedogonata 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ8—06 

Arundinaria faberi 

Hangzhou Bamboo Garden 

ZHUGE Qiang 

HZK—5 

Brachystachyum densiflorum 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

Nil 3—04 

Phyllostachys edulis (outgroup) 

Bamboo Garden of Nanjing Forestry University 

DING Yu-Long 

NJ21—03 


Note: ©The vouchers were stored at the herbarium of Nanjing Forestry University 


Results and discussion 

The length of the entire ITS region (including 5.8S) in all taxa 
sampled ranged from 596 to 623 bp, with ITS1 ranging from 215 
to 225 bp, 5.8S 163 bp for all species, and ITS2 ranging from 
213 to 231 bp, and the (G+C) contents in ITS2 (ranged from 
72.3% to 76.9%) are slightly higher than those in ITS 1 (ranged 
from 69.6% to 74.5%). In the data matrix, the levels of 
potentially informative nucleotide site (14.3%) of ITS1 are 
somewhat higher than those of ITS2 (6.7%). The similar pattern 
was also found in Krigia (Kim & Jansen 1994), Populus (Shi et 
al. 2001). The aligned tmL-F data matrix has 116 variable 
characters and 15 (1.5%) potentially informative characters out 
of 973 characters. The sequences are, typically for chloroplast 
DNA, rich in AT (66.6%). Of the 1 658 aligned ITS and trnL-F 
combined sequences, 1405 (84.7%) are constant, 186 (11.2%) 
uninformative and 67 (4.1%) are informative (Table 2). 

Heuristic searches including only ITS characters generated 
three most parsimonious trees requiring 205 evolutionary steps 
(strict consensus in Fig.l) with a consistency index (Cl) (0.8293) 
and retention index (RI) (0.7083). The ITS data indicates that the 
species sampled in Arundinaria and the related genera are 
monophyletic and divided into two clades with stronger support 
(bootstrap values >50%). One clade included A. oleosa, A. 
hsienchuensis, A. chino, A. amara, A. yixingensis, A. amabilis, A. 
fortunei and A. pygmaea, and the other contained A. gramineus, 
A. fargesii, A. faberi, A. hupehense, Pseudosasa japonica cv. 
Tsutsumiana, P. japonica, Brachystachyum densiflorum, A. 
oedogonata and A. sulcata. Therefore, it is likely that all these 
taxa sampled belong to Arundinaria , which should include the 
related genera ( Pleioblastus , Pseudosasa, Oligostachyum. 
Bashania, Clavinodum). The result also showes that there is a 
close relationship between A. gramineus and A. fargesii, 


Pseudosasa japonica cv. Tsutsumiana and P. japonica, A. sulcata, 
Brachystachyum densiflorum and A. oedogonata (99%, 100% 
and 87% bootstrap support, respectively) . Furthermore, it is 
shown from the phylogenetic tree that A. pygmaea and A. 
fortunei have a close relationship, and are a sister branch to 
Pleioblastus bamboos. 


Table 2. ITS region and trnL-F intergenic spacer (IGS) sequences 
characteristics of the bamboos sampled 


Parameter 

ITS 

trnL-F 

combined 

Number of taxa 

18 


18 

18 

Sequence length range (bp) 

596-623 

868-963 

1467-1560 

Total aligned matrix length (bp) 

669 


973 

1658 

G+C content (%) 

69.7 


33.4 

48.1 

Number of constant sites, % 

519(77.6%) 

842 

(86.5%) 

1405 (84.7%) 

Number of variable sites , % 

96(14.3%) 

116 

(11.9%) 

186(11.2%) 

Number of uninformative sites , % 

42 (91.9%) 

101 

(98.5%) 

119(95.9%) 

Number of informative sites , % 

54 (8.1%) 

15 

(1.5%) 

67(4.1%) 


The phylogenetic analysis by trnL-F IGS region sequences 
produce one most parsimonious tree (Fig. 2) of 139 steps length, 
with consistency index (Cl) (0.9495) and retention index (RI) 
(0.7500). The topology of the tree of the trnL-F dataset is similar 
to the ITS tree, and also showes a monophyly of the species 
sampled in Arundinaria and the related genera which are fell into 
two clades. One clade includes the same species with similar 
relationships as in the ITS tree, except that the position of A. 
amara and A. yixingensis are replaced by A. gramineus, A. 
sulcata, P. japonica cv. Tsutsumiana and P. japonica. 

The maximum parsimony analysises of the combined data of 
ITS region and trnL-F IGS region generate one most 
parsimonious tree (Fig. 3) of 332 steps length, with a consistency 
index (Cl) (0.8404), a retention index (RI) (0.6936) and rescaled 
consistency index (RC) (0.5829). The phylogenetic tree from 
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both sets of data is broadly similar in the topologies of the trees 
from the individual data, showing a monophyly of the species 
sampled in Arundinaria and the related genera (including 
Pleioblastus, Pseudosasa , Oligostachyum, Bashania, 
Clavinodum, etc.) which are divided into two clades. One clade 


comprises Pseudosasa japonica cv.Tsutsumiana, P. japonica, A. 
yixingensis, A. amara, A. oleosa, A. hsienchuensis, A. amabilis, A. 
chino, A. fortunei and A. pygmaea, and the other comprises 
Brachystachyum densiflorum, A.oedogonata, A. sulcata, A. 
hupehense, A. fargesii, A. gramineus and A. faberi. 



A. oleosa 
A. hsienchuensis 
A. amabilis 
A. chino 
A. amara 
A, yixingensis 
A. pygmaea 
A. fortunei 
A. gramineus 
A. fargesii 
A. faberi 

A. hupehense 

P. japonica cv. Tsutsumiana 
P. japonica 

B. densiflorum 
A oedogonata 
A. sulcata 

Ph. edulis 


Fig. 1 Strict consensus tree of the 3 most parsimonious trees from the entire ITS sequences with gaps coded as missing. 

Tree length, 205 steps; Cl, 0.8293; RI, 0.7083. Numbers above lines represent bootstrap values in 1 000 replicates. 



A. gramineus 

P. japonica cv. Tsutsumiana 

P. japonica 

A. oleosa 

A. hsienchuensis 

A. amabilis 

A. sulcata 

A. pygmaea 

A. chino 

A. fortunei 

A. hupehense 

A. yixingensis 

A. fargesii 

A. faberi 

B. densiflorum 
A. oedogonata 
A. amara 

Ph. edulis 


Fig. 2 The most parsimonious tree from the trnL-F IGS region sequences with gaps coded as missing. 

Tree length, 139 steps; Cl, 0.9495; RI, 0.7500. Numbers above lines represent bootstrap values in 1 000 replicates. 
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P. japonica cv. Tsulsumiana 
P. japonica 
A. yixingensis 
A. amara 
A. oleosa 
A. hsienchuensis 
A. amabilis 
A. pygmaea 
A. chino 

A. fortunei 

B. densiflorum 
A. oedogonata 
A. sulcata 

A. hupehense 
A. fargesii 
A. gramineus 
A. faberi 
Ph. edulis 


Fig. 3 The most parsimonious tree from the combined data of ITS region and trnL-F IGS region sequences with gaps coded as missing. 

Tree length, 332 steps; Cl, 0.8404; RI, 0.6936. Numbers above lines represent bootstrap values in 1 000 replicates. 


Our study indicated that the two regions of ITS and trnL-F 
IGS provide an important clue for the establishment of 
phylogenetic relationships of the bamboos. The ITS sequence is 
not only alignable at intergeneric level in Arundinaria, but it is 
also phylogenetically useful at this level. Comparison with 
trnL-F IGS sequence, the ITS region provides the higher number 
of parsimony informative characters, and the interspecific 
variation of the ITS sequence is higher than that of the trnL-F 
IGS sequence. The topologies of the trees generate from 
individual data and the combined data are similar. The results are 
significant in confirming aspects of previous morphological 
analyses (Yang & Zhao 1993). 
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